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ON THE DYNAMIC BEHAVIOUR OF A ROTAMETER 


J. F. Ury 
Faculty of Mechanical Engineering, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The position of the float in a variable-area flowmeter (rotameter) for liquids is determined 
by the velocity of the liquid in the annulus around its measuring edge. Pulsations in this 
velocity bring about periodic variations of the float position. The equilibrium of the forces 
exerted on the float by the relative velocity of the pulsating liquid with the forces of gravity 
and inertia of the float is expressed by a non-linear differential equation, of a form not 
suitable for practical application. For small variations, however, a linear second-order 
approximation is obtained, readily predicting the performance of the meter with regard 
to magnitude and phase shift of the instantaneous displacement. As long as this approxi- 


mation holds, no error of mean measured value should exist. Test results fully confirm 
these predictions. 


In this paper an equation is derived for pulsating liquid flow through a variable 
area flowmeter (rotameter). It is obtained by analogy to the well known rotameter 
equation for steady flow. But it is a non-linear differential equation, and the flow 
coefficients appearing in it cannot be determined by direct calibration, as is pos- 
sible for steady flow. Numerical solutions can, therefore, be found only in special 
cases in which it seems possible to use linear approximations, together with esti- 
mated values of the coefficients. 


This paper deals only with one of them: the case of small accelerations. For these, 
performance calculations were carried out and compared with the results of tests, 
concerning two main problems: 


a) One of foremost theoretical interest: the instantaneous response of the meter 
to sinusoidal forcing, expressed in terms of gain and phase shift. 


b) One of higher importance for the practical use of the meter: the rotameter er- 
ror, i.e. the difference between mean indication and true average flow. 


THE FLOWMETER EQUATION 


An exact analytical treatment of rotameter performance would have to be based on 
the fundamental equations of hydrodynamics, integrated over all surfaces and vol- 
umes concerned. Interesting steps in this direction were taken by Kramers}; but it 
seems impossible, at least at present, to pursue this course to the end by considering 
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to full extent all the various effects actually present, including i.a. the varying dis- 
placement of the flow restriction itself, unequal distribution of velocities, viscous 
forces etc. In these circumstances, empirical correction factors must in any case be 
applied, at a certain stage, to the results of theoretical reasoning. Some doubt may 
therefore arise whether the extra work spent in deriving and using such (necessarily 
rather complicated) forms of equations would be justified. 

In the following it was attempted, instead, to obtain an equation for pulsating 
liquid flow by analogy to the steady flow equation derived by Head* on the basis of 
a highly simplified theory. Theoretically such an approach may appear much in- 
ferior, and certainly an equation arrived at in this way should not be expected to 
give information on what is actually occurring in the flow around the float at each 
particular instant. But it will be shown that it can be used, in special cases at least, 
to predict the overall performance of the meter witha satisfactory degree of accuracy. 


The equation for steady flow is4: 
We = K(a? — 1) a /F Dy Fup (1) 


where K is an empirical flow coefficient taking into account, besides the contraction 
of the jet, all influences not explicitly considered in the simplified theory, e.g. vis- 
cous drag, unequal distribution of velocities etc. 

With w,, expressed by the water velocity x,,; in the actual smallest section of the 
jet around the measuring edge 


Wst a PajXwj = parC (a? = OFC os 


eq. (1) can be written in the form: 
st = 0 (la) 


stating the equilibrium of forces acting on the float. The first term is the force due to, 
and acting in the direction of, the water velocity x,,; (upwards), and the second — 
the restoring force, i.e. the apparent weight of the float in water (acting downwards). 

This fundamental equation of forces must also hold for pulsating flow, provided 
values corresponding to dynamic conditions are used, for the relative velocity of 
the water across the measuring edge: 


_ _ W payx 


X= papCofe ergy —x (7a) (see Appendix 1) 


and for the restoring force: 


Fan ite E + ier: ;| (8a) (see Appendix 2) 
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In addition, the fact that the velocity force on the float always acts in the same 
direction as the velocity causing it, must be expressed mathematically by the term: 


SAR 


the flow coefficient K being given the same sign as x,,. 


The force equation then assumes the form: 


a gti) Py amy 2) 


showing that the restoring force now may become zero and even negative, if acceler- 
ation amplitudes are sufficiently large. It changes its sign simultaneously with x,; 


and, therefore, with K. 
ep OOD Sapeeea ay ' 
oh C.p. ay J |KI Fayn p (2a) 


with the term under the radical essentially positive. 


It follows that: 


By equating with eq. (7a) and rearranging, the flowmeter equation for pulsating 
flow is obtained: 


aa 
w = Ka? — 1a i Dy J Fin 2 + pagk(l + C.(@2- 1D) @) 


Me 
or, for a conical tube with Dy = D, and hence « = —— =1 + 6x: 


w = 2Kdx [1 + sox (1+ 8s) o,f Rag + 
+ payx [1 + 2C,dx(1 + oe) (3a) 


of which eq. (1) or ‘ 
iw 2 Kox{1 # 5ox( yf SS) v2" D, Fa (1b) 


is the particular case for steady flow. 


For practical purposes, this equation must be further simplified. The rotameter 
scale is determined by calibration under steady flow. It is found to be nearly linear. 
Thus little error will be introduced by assuming that the flow coefficient K varies so 
that 


K 
K [1 " 5 5x(1 r 2) re iE — const. (3b) 
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In an analogous way it may be assumed that 


2C 0 (1 + =) = K, = Const. (3c) 


so that the rotameter equation can be written in the form: 


K ; : 
w= stk RY Fayap ot pa,K>x]| + pa sx (4) 
1 


It should be borne in mind that eq. (1) was derived with the aid of Bernoulli's 
theorem for steady flow. For pulsating flow, an acceleration term appears, containing 


the integral es dx taken over the axial length of the flow restriction. This also 


has to be taken into consideration in the flow coefficient K,, in addition to the 
other effects included in K for steady flow. 


For pulsation amplitudes exceeding in magnitude the average water velocity, the 
flow relative to the float will be directed downwards during part of the cycle. 


If, then, K,, (>0) and K,, (<0) are introduced to denote the respective average 
time-values of K, during periods of upward and downward relative flow, and using 
the ratio: 

2 hBaake 

[Kaul ” 


Ri 


and eq. (8a), the equation assumes the final form: 


ee ee Peer 
w=x!K,./F thie So someill tdy Br ee ates 
where only real values of the roots have to be considered. 


To be exact, the coefficient K», too, should be assumed to depend on the direction 
of relative flow. But the term in which it appears is relatively unimportant. Therefore 
it was considered that the assumption of a mean value for both directions will not 
introduce an appreciable error in the results. 


Tests indicated that, at any given frequency, the forcing amplitude can be varied 
within rather wide limits, with no significant effect on the gain of the system. This 
leads to the conclusion that the coefficients K; and K> are, at least at first approxi- 
mation, independent of the forcing amplitude, whereas they will vary with frequency®. 


The solution of the equation depends on the term: 
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which has to be treated in different ways depending on the magnitude of the acceler- 
ation term under the radicals. 


In the present paper only the case of small accelerations is dealt with. 


APPROXIMATION FOR SMALL ACCELERATIONS 


It can be easily shown that the acceleration term can be approximated by a Fourier 
Series with cosine terms only: 


Pr 
—/_—=p,coswt + Y p,cos (not + &,) 
Prp-pg *' » 


and that the coefficients of the higher harmonics are small compared with py. 


A 


This means that as long as p, = Ps % co? ~ is small compared with 1, the term 
"OE a 
ahs San (er aa) 


becomes approximately, to the first order of small variations, 


adraa ea il it 
ptt ena + ppp 29 
and, therefore: 


ays ae Ps K2x| - : 5 
W = x{ Ky veeAll + ppaaip *) + pi 2%] + pax (5) 
which, by substituting x = x,, + Ax and accordingly w,, = %,K1./Fp, assumes 
the form: 

Ky 


Stay 5 VF ab (Ax) + pa,(Kz + -. (Ax) + At LIF ap (Ax) +, 


W — Wm = Xn {a 


5 (5! er sieef) (Ax) +- pa,K(A%)| + 


m 


= x,,{S(AX) + T (Ax) + U(Ax) + (Ax) [V(AX) + W(Ax)]} (6) 


where the terms: 


were introduced for brevity. 
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The last two terms of the equation, containing products of variations supposed 


to be small, may be neglected with small error. 


The equation then becomes approximately: 
Wo — Wy & Xm{S(Ax) + T (Ax) + U(Ax)} (6a) 


which is a linear second-order equation. As long as K, and K, can be assumed ap- 
proximately constant, this equation has constant coefficients and defines a physical 
system of second order2, having the constants: 


aT: P45 (K2Xm + 1) 
eh’ SAFE 


Characteristic time: T= 


Undamped natural frequency: f, = <—,/—- = 


1 ju sa behcara 
Xm PF 


T KyXnt+1) [2% vy 
Damping ratio: ¢ = 2tf, = + gs = As ( ) Gg p ee ; Pp) 
Liat Fy, 3 - 


by which the response of the system with regard to amplitude and phase shift is 
determined. 


As stated above, the most essential point concerning the practical use of a rota- 
meter is, however, the determination of the error of the mean indication, i.e. the 
deviation of the apparent mean flow from the true average flow, or Wy — Way 


From Ma Wait (w i Wav) FF. (Win a Wav) 
follows the time‘average: w —w,, = (W— W,,) — (Wn — Wap) 
And since by definition w—w,, =0, it follows that 


Wn 7 Wav = Ky J Fup Xe = Xa) 1 a (w— Wr): 


Since the time average of all terms on the right side of eq. (6a) is zero, it can be con- 
cluded that the measuring error is negligibly small as long as this approximation 
holds, i.e. as long as variations are small. This is in agreement with Heads. 


MEASUREMENTS 


Measurements for verification of the theory were carried out at the Instrumentation 
Laboratory, Department of Technical Physics, Technological University, Delft 
(Holland). The test installation used was devised and built by Schneiders and Kramers, 
who were also the first to use it for measurements?.3, 
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Figure 1 


It is shown schematically in Figure 1. Water was used as the flowing medium. 
Its level is held constant at 1, and the flow adjusted by cock 2, and measured by 
rotameter 3. The average flow may also be measured by a measuring flask and 
stop watch at discharge 10. 

Pulsations are generated by metal bellows 4, driven by a d.c. electric motor with 


the aid of crank mechanism 5, the radius of which is adjustable. The indication of 
the rotameter is compared with that of a magnetic flowmeter 6, which can be as- 
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sumed to measure true instantaneous velocities!. It comprises a field{coil with ex- 
citation by alternating current of about 400 cps, electrodes, amplifiers, and filters, 
and one beam of a double-beam cathode ray oscillograph 7. The other beam is deflec- 
ted by the output voltage of photocell 8, which drops when the light beam 9 is inter- 
cepted by the rising float. In this way the phase shift between water velocity ard 
float movement is indicated on the screen. Figures 4, 5 and 6 show examples of oscil- 
lograms taken at different frequencies. 


Figure 4 Figure 5 Figure 6 
f=3.1 cps f=F2- cps f = 17.4 cps 


In the summer months of 1957 the author was given the privilege of working as 
a guest in that laboratory and carrying out a number of tests with this installation, in 
order to verify the equations with regard to correctness and limitations. Only test 
data obtained during that period by the author himself are included in this paper. 


Since preliminary tests had indicated that change of crank radius had no signi- 
ficant influence on response, all final tests were carried out with a constant crank 
radius of 1.25 mm, corresponding to 3.86 cm3 full-stroke water displacement. Several 
series of measurements were carried out, each of them at a constant opening of cock 
2 corresponding to a certain constant steady flow. One of them, at 45% of maximum 
rotameter flow, is given here as a typical example. 


Figure 2 shows how water velocity amplitudes vary with frequency. In the range 
of measurements, two resonance peaks exist. It is believed that this is due to the low 
natural frequency of the elastic bellows used as pulsator. Nearly identical response 
is obtained when the float is removed, showing that the response of the partial system 
pulsator-water is not altered significantly by the insertion of the float. In Figure 3 
the deviation of actual average flow from steady flow is shown as function of frequency. 
In this respect the differences due to the presence of the float are obvious, concerning 
the magnitude of minimum throughput as well as the frequency at which it 
occurs. 
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BEHAVIOUR OF THE ROTAMETER 
(a) Response to sinusoidal forcing: 

The rotameter used was Fischer and Porter, type B4-21-10/70, with a stainless steel 
floattype BSVT 44. Figure 7 shows amplitudes of float movement xX, expressed in percent 
of maximum rotameter flow. This meter’ is graduated in percent of maximum flow, 
100% corresponding to 187 1/h water flow. The zero mark is not indicated, but the 
average increment of tube length per scale division corresponds to x;99 = 0.256 m, 
equivalent total tube length between 0 and 100%. Further data of the meter are: 


Do = Dp=13:%107 4m, agi a= 1,33,x 105 4m? 


. — Pr00 = Do = (0.785 m7} 


D, X 100 


pr = 8000 kg m™*, p=1000kgm°, F =0.126 Newton. 
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Thus for small variations the characteristic constants of the system are as fol- 
lows: 
0.0119 
t= 
Ky, 


(KiXm +1) (sec) 


_ 0.0247 Ky% m+! 
Kyu sh Be 


¢ 


which, besides being functions of x,,, are also dependent on the flow coefficients 
K, and K, which vary with frequency®. 


If it is assumed, however, that values of K and C, valid for steady flow may still 
be used approximately for small variations, K, = 0.0178 and K, = 1.04 are ob- 
tained when mean values (acc. to Head)4 of K = 0.6 and C, = 0.63, together with 
x = 0.128 m (one-half of total length of scale), are substituted in eqs. (3b) and (3c). 


This results in: 
1 = 0.668 (1+ 1.04x,,) and ¢ = 1.39 oP 


wherefrom Table I was calculated: 


TABLE I 
Average flow mz t r a 
Max. rot. flow 
(%) . (m) (sec) (cps) 
0 0 0.688 CO ee) 
45 0.115 0.742 1.944 4,53 
100 0.256 0.846 1.306 3.47 


The corresponding curves for gain and phase shift are plotted in Figure 8, together 
with the results of the measurements. The curves define a band of response of the 
rotameter in which the test points for every mean flow should lic. The lower limit, 
the first-order system with t= 0.668 sec, is of course not attainable in practice, 
since then the lower stop would interfere with free oscillation of the float. 


It is seen that agreement between the approximate equations and the tests iS 
fairly close not only at frequencies below 1.5 cps, where p; is small (see Figure 7), 
but even at frequencies between 3 and 5 cps, where pj > 1. That all test points © 
around 2 cps lie distinctly below the approximate curve would agree with the theo- 
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retical postulate that K, should decrease®, and, therefore, t increase with frequency. 
That the agreement is better again at somewhat higher frequencies, and remains so 
up to about 6 cps, may be due to a combination of mutually counteractive effects, 
and cannot be explained by the simplified theory. 

At still higher frequencies no agreement is found, although in that case p, dec- 
reases and the approximation should hold again. Two main causes appear to be 
responsible for this contradiction: 


1. The distortion of the sinusoidal forcing by higher harmonics in the frequ- 
ency region around 8 cps, which is recognisable, for example, in Figure 5. 

2. The limited accuracy with which the very small amplitudes of float movement 
occurring at still higher frequencies can be measured. 


(6) The rotameter error: 
Throughout the whole frequency range in which p, was found small (say smaller 
than 0.4, corresponding, in the test installation used, to f < 2.5 cps and f > 9.5 
cps approximately; see Figure 7 as an example), no significant rotameter error 
(i.e. one exceeding 2% of the full scale) was measured, in complete agreement with 
the approximate equation for small variations. 


CONCLUSION 


(a) The non-linear differential equation for pulsating liquid flow through a rota- 
meter can be approximated, for small variations, by a linear one of second order 
with constant coefficients. Tests showed that the response of the meter with regard 
to gain and phase shift is represented fairly closely by that same approximate equa- 
tion throughout the whole range of frequencies in which purely sinusoidal forcing 
can be assumed to exist, even for frequencies at which there is no theoretical jus- 
tification for this agreement. 

(b) As long as variations are small, justifying the use of the linear approximation, 
the error of mean float position does not exceed the limits of accuracy of the meter 
for steady flow, estimated as + 2% of the full scale. 


APPENDIX | 
The relative velocity x,; of the water through the actual smallest section across 


the measuring edge is found by superposition of two velocities: 
a) x;,, due to the motion of the float in water assuined to be at rest: 


af } 1 
xy =—X E OF (ha LA), i) | 


b) x, due to the pulsating flow of the water along the float held in place: 
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Thus, by superposition: 


Swen 


g op, aL, wo x 
Np Mey eae py (7), 
IAT EY TO EY, FON) 
or,, with pa x, =W, 
W — payx 


(7a), 


a ~ parC, (a2 = A)usi 


which can be positive, zero, or negative, depending on instantaneous values of w 
and x. 


APPENDIX 2 


The restoring force for pulsating flow was calculated as follows: 


Let the volume of the float be V,, its density p,. The water streamlines passing the 
float cannot follow its contours exactly. Some deadwater zones are formed, so that 
a total volume Vip > V; (represented by the shaded area of Figure 9) is accelerated 
by the pulsating flow, Assume, at first, the float to have been removed, but the 
shape of this “body” of volume V; retained by an infinitesimally thin, but rigid, 
shell of negligible mass, entirely filled with water. Then the external force required 
to generate in it the acceleration x would be: 


Now assume the float to be inserted. Then, since Vij is constant, the same volume 
V, of water will be driven out. If such an insertion is made in a region with a field of 


Aypothetic shell 


Figure 9 
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acceleration g, the external force required to balance the float is, according to hyd- 
rostatics: 
Vs (py — p)g 


But here, in this “body”, the total acceleration is g + x instead of g. Therefore the 
additional balancing force is proportionately increased to: 


+X x 
F, =Vs(0g — 99 “= =Ve(,—p)g(1+-) 


The total external force which has to be exerted in this case on the moving “body” 
by the pulsating flow of water passing it is, therefore, 


Fyn = Fy + Fy = Vz pX + Ve (py, — p)(g + X) 


Vesa Pls 


= V;(py — p)g + Vjp, [1+ Woy ig ? 
f f 


Denoting the volume of deadwater accelerated together with the float, per unit 
volume of float, by 


V/=V 
elarreae 
; f 
we obtain 
a Ps P\xX 
Fan = Fu [1+ (14+ 97) (8) 


At present nothing is known about the values which the ratio m may assume under 
different conditions of pulsating flow, and tests on this point are under way. It 
seems, however, permissible to assume that as long as a stainless steel float is used, 
with 


the second term in parentheses will be small compared with 1, and may therefore 
be neglected at first approximation, This gives: 


Fan © Fa[ 1 + ae (8a) 


This form was used in the derivation of the rotameter equation in this paper, 
which deals exclusively with the performance of a stainless steel float. 

Some modification may prove necessary in the case of floats of lower densities 
(aluminium, plastics etc.). It is intended to deal with these problems at a later date 
in a separate paper. 
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NOTATION 
a m2 section 
a, a >»? of meter tube at zero of scale 
A100 * *»» *»* of meter tube at full scale 
a, re area of annulus around measuring edge of 
float 
ay “ largest section of float 
a; a = C,a,, actual smallest area of water flow 
a, 4 = a2a,, section of meter tube at height x 
C, nondim. coefficient of contraction 
D om diameter 
yp: :? ** °” of meter tube at zero of scale 
Dio0 7 ae *» at full scale 
D; 3 largest diameter of float 
D . tube diameter at height of measuring edge 
(m or rotameter error 
E jy of max. 
\rot. flow 
et Newton _ restoring force for steady flow 
oe et Ey 7 *” for. pulsating flow 
cps frequency 


Se Pm undamped natural frequency 
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K, K,, K, 


Wav 


Ps 


nondim. 


99 
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rot. flow 


99 
99 


9° 


m ss" 
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nondim. 
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flow coefficients 

flow coefficient for upward relative flow 
ss for downward relative flow 

Fourier coefficients 

parameter as defined in text 

ratio of flow coefficients 


factors of terms in rotameter equation as 
defined in text 


mass flow 


o> 3° 


, average 


” , corresponding to steady flow with 
float at x,, 


position of float 


position of float corresponding to w,, 
mean position of pulsating float 


amplitude of sinusoidal variation of float 
position 


velocity of float 

velocity of water 

amplitude of sinusoidal variation of velocity 
acceleration of float 

acceleration of water 


= D,/D,, ratio of tube diameters 


shape coefficient of meter tube 
phase shift angle 
variation about mean value 
damping ratio 
characteristic time 
density of water 

of float 


angular velocity 
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DETECTION AND DETERMINATION OF KEROSENE AND GAS 
OIL CONTAMINATION OF M.T. GASOLINE 
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ABSTRACT 


The effect of small admixtures (up to 2%) of kerosene and gas oi! on the properties of 
local M.T. gasoline was studied. Three basic methods were used: distillation, bromine 
number and paper chromatography. Semimicrodistillation of the Engler residue (6%) 
accentuates considerably the increase in F.B.P. The bromine number method was found 
unsuitable owing to the mutually counteractive effects of the decrease in olefin content 
and increase in sulphur content. Paper chromatography permits both identification and 
quantitative determination of the contaminant. 


INTRODUCTION 


The problem of detection of small admixtures of kerosene and gas oil in M.T. gasoline 
is still unsolved despite its recently acquired world-wide importance; in fact, little 
material is available on this subject, with the single exception of a chromatographic 
method for detection and determination of contamination of aviation gasoline!. 

Local M.T. gasoline is an 83-octane product, defined according to I.S. 90. Small 
admixtures (up to 2°%) of kerosene or gas oil affect some of its properties, but the 
fluctuation is too narrow to justify rejection in accordance with I.S. 90. 

The present study was undertaken with a view to developing analytical methods 
permitting the determination of kerosene and gas oil in the gasoline in amounts 
below 2%. 
GENERAL 
Utilisation of the conventional fuel testing methods for determining small admix- 
tures of kerosene in gasoline had been studied at this laboratory. These methods 
require a standard for comparison and can be used for amounts of over 3% kero- 
sene. Methods based on adulteration of kerosene, with the adulterants serving as 
tracers, had been studied both at the Shell Laboratories in England? and this la- 
boratory. The most common adulterants used are: 

1. Furfural — permits detection of 2% kerosene in gasoline, but the colouring 
substance and inhibitors interfere with the analysis. 

2. Red detector F.L. — unsuitable for the detection of small amounts of kero- 
sene; its colour is noticeable even at low concentrations. 
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3. Fluorescein — detectable in small amounts, but its solubility in kerosene is 
too low to give a stable solution for more than a few days. 


4. Radioactive tracers — detectable even in the minutest amounts, but fairly 
expensive. In addition, the possibility of dangerous concentration levels in the 
containers should be borne in mind, so that their suitability cannot be ascertained 
without further research. A study on this subject? was published a few years ago, 
but the authors were unable to obtain it. 


In the present study, three basic methods were used, namely: 


1. Distillation. 
Bromine number. 
3. Chromatography. 


1. Distillation 


Since the refining of crude oil is based on various distillation processes, it could 
be assumed that these methods can be utilised for the identification of heavier 
hydrocarbons. In fact, there are numerous methods based on their separation by 
distillation**© but these cannot yield satisfactory results in view of the fairly wide 
overlapping, in the local products, between the heavier fractions of gasoline and 
lighter fractions of kerosene. Table I shows that about 50% of the kerosene hydro- 
carbons distil over up to 188°C, i.e. within the distillation range of gasoline. It is 
obvious that in these circumstances a 2° admixture of kerosene could not be de- 
termined by the ordinary Engler distillation test prescribed by A.S.T.M. Standard 
D-867. 


An attempt was made accordingly to obtain better separation of the hydrocar- 


bons by the following methods: 
a) Distillation according to A.S.T.M. Standard D-2857. 


b) Preliminary distillation up to 80% in a column charged with Raschig rings, 
and treatment of the remaining 20% by the D-86 and D-285 methods. 


c) Preliminary distillation of the gasoline by the D-86 method up to 94%, 
concentrating the greater part of the contaminants in the residue, and treating the 
remaining 6% by a semimicro method. 


2. Bromine number 


Since locally marketed gasoline is produced by a cracking process, its bromine 
number is high (above 30) while the kerosene, produced by straight-run distillation, 
has a bromine number near zero. The reduction in bromine number due to’ 2% 
kerosene could be too small to be detected by the ordinary method, so the test was 
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TABLE I 


Comparison of distillation ranges of local gasoline and kerosene 


Gasoline Kerosene 

WBE. 43°C # 1856°C 
10% 60 164 
20% 69 170 
30% 79 176 
40% 89 181 
50% 99, "ta. 188 
60% 108 | 196 
10% 117 206 
80% 128 218 
907, 146 231 
eBaE: 185 244 


carried out on the residue. In this case, certain factors likely to affect the bromine 
number of the residue have to be taken into consideration, such as the possibility 
of cracking during distillation and the concentration of sulphur compounds likely 
to react with the bromine, such as:8 


R»S + Br2 - RS. Bro 
R,S. Br, + H,0 > R,SO + 2HBr 


R2»S + Bro > (R,S*. Br)Br- 
(R2S*. Br) Br~ + H,O > R2gO + 2HBr 


3. Paper chromatography 


Krynitsky and Garrett! developed a chromatographic method for the detection 
of aviation kerosene in aviation gasoline products with non-overlapping  distil- 
lation ranges. The method is based on immersing a strip of filter paper in the tested 
gasoline and subjecting it to a current of air, causing the absorbed liquid to evapo- 
rate. The lighter fractions evaporate more rapidly, so that the gasoline contami- 
nated with heavier compounds shows a higher level of diffusion than pure gasoline. 
The method developed by the authors, in an attempt to adapt this technique for 
motor gasoline with a high degree of overlapping (see Table I), permits quantitative 
determination of admixtures in amounts exceeding 0.5°% as well as identification 
of the contaminant (kerosene or gas oil). 
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Figure 1 
F.B.P. vs. kerosene content (semimicrodistillation method) 


EXPERIMENTAL AND RESULTS 


A. Distillation 


Comparisons of the D-86 and D-285 methods and semimicrodistillation were car- 
ried out. In the latter case, the procedure was as follows: 


Preliminary distillation was carried out by the D-86 method up to 94%. 5 ml of 
the residue were treated in a 20 ml. Claisen flask and 10 ml. Wurtz flask. Figure 1 
shows the increases in final boiling point as functions of the kerosene content. 


The semimicro method accentuates the difference between the final boiling points, 
with the most pronounced effect obtained by using a Wurtz flask—a 20° increase 
for 2% kerosene and 35° for 4% kerosene. 
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F.B.P. vs, contaminant content (D-86 method) 
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Figure 3 
Distillation residue vs. contaminant content (D-86 method) 
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TABLE II 


Bromine number vs. sulphur content 
a ce a ie eS ae I ee eee ciel i 


A.S.T.M. D-86 Distillation Vacuum Distillation 
Sulphur Bromine Sulphur Bromine 
(Y) number (%) number 
Frac- Pure Gasoline Pure Gasoline Pure Gasoline Pure Gasoline 
tion gasoline + 2% gasoline + 2% |gasoline + 2% — gasoline +.2% 
(ES) kerosene kerosene kerosene kerosene 
0—20 63.6 63.5 
20—40 45.7 43.9 
40—60 32.4 30.9 
60—80 0.14 24.7 2G 
80—93 0.20 0.19 26.0 25.4 
Residue 0.38 0.28 37.9 33.0 0.67 0.71 42.8 36.4 


The residue obtained for a 1% admixture of gas oil has a final boiling point of 
275°C in the Wurtz flask, i.e. an increase of 72°. 


In another test, distillation by the D-285 method was carried out up to the final 
boiling point and the residue treated in a Wurtz flask. For pure gasoline the residue 
showed a boiling point of 223°C; with 2°% kerosene — 230°C. The small difference 
is explained by the heavy reflux, resulting in the concentration of heavy fractions 
in the residue. 


Conclusions 

1. The D-86 method is the most effective means of accentuating the presence 
of foreign substances affecting the final boiling point (Figure 2). 

2. This method indicates the presence of gas oil by the greater amount of re- 
sidue compared with kerosene (Figure 3). 


3. Preliminary distillation up to 94% by the D-86 method and semimicrodis- 
tillation in a Wurtz flask give a pronounced increase in final boiling point which, 
even for 4% kerosene, is fairly close to that of the pure kerosene. An admixture 
of 1% gas oil gives an extremely high boiling point. 


B. Bromine Number 


Ordinary (D-86) and vacuum distillation residues of M.T. gasoline containing 
1, 2, 4 and 6% kerosene were tested for bromine number. Results were of the same 
order in both series, the conclusion being that ordinary distillation does not pro- 
duce cracking. In addition, the sulphur content and bromine number were deter- 
mined for the residue and various fractions by the I.P. 61 method.9 Table II shows 
that the higher the fraction, the lower the bromine number which increases again 
for the residue owing to the increase in sulphur content. Considering the mutually 
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counteractive effects of reduced olefin content and increased sulphur content, this 
method is unsuitable for the detection of kerosene. (In addition, its reproducibi- 
lity according to I.P. 129 is 2). 


C. Paper Chromatography 
a) Equipment and reagents. 
1. Ventilated chamber (Figure 4). 
2. Beaker, tall type, 100 ml, 45 mm diameter, 75 mm height (Figure 5). 
3. Cover—aluminium foil, slotted in the middle. 
4 


Strip of Whatman No. | chromatographic paper, 205 mm length. Zero line 
at the slot (35 mm above lower end, depth of immersion 10 mm). 


5. Dye—5% solution of orange dye in benzene. 


b) Procedure. 


Ten drops of the dye solution are added to 80 ml of the tested gasoline in the 
beaker. The paper strip is suspended over it and immersed in the gasoline through 
the slot of the cover. Air, at the rate of 30 ft/min, is blown through the chamber 
for 2 hours. The strip is then withdrawn and dried for 15 minutes in the same air 
eurrent, after which the level of diffusion is measured. 


The method was arrived at as a result of a series of experiments with a view to 
establishing optimum conditions for the following details: (1) dye, (2) shape of 
paper strip, (3) rate and duration of ventilation, (4) effect of ambient tempera- 
ture and humidity. The two last-named factors have considerable influence on 
the diffusion level. Table III gives an example of a test showing that the level of 
diffusion rises with humidity owing to the reduced rate of evaporation. 


TABLE It 


Effect of ambient temperature and humidity 


Humidity (°%) Temperature (°C) Level of diffusion (mm) 


56 24 33 
64 22 39 


144 


M. BABITZ,S. BENTUR AND A. ROCKER 


Figure 4 
Ventilated chamber 


Figure 5 
Beaker with immersed paper strip 
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c) Tests 
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Figure 6 


Effect of gasoline properties on level of diffusion 


1. Effect of gasoline properties on level of diffusion 


Samples were varied according to specific gravity and final boiling point. Results 
are in Table IV, showing that: 
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i) The heavier the fraction, the higher the level of diffusion (Samples 1, 2 
and 5, Figure 6). 


ii) Similar samples (No. 2, 3 and 4) gave similar levels of diffusion. 


TABLE IV 


Levels of diffusion of various gasolines 


Sample No. Specific gravity sae Level aa 
1 (straight run) 0.7032 145 11 
2 0.7312 174 31 
3 0.7318 175 33 
4 0.7317 178 33 
5 0.7325 188 46 
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Figure 7 : 
Effect of kerosene content on level of diffusion 
(1) Pure gasoline. (2) 0.5% kerosene. (3) 1% kerosene. (4) 2% kerosene. 
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_ Figure 8 
Level of diffusion vs. kerosene content 


2. Effect of kerosene content on level of diffusion. 

Results are given in Figures 7 and 8, showing that the level of diffusion rises 
with increasing kerosene content. Pure gasoline showed in all cases a high con- 
centration of the dye at the upper edge of the chromatogram. The presence of 
kerosene is indicated not only by the higher level of diffusion, but also by the even 
distribution of the dye over the chromatogram, resulting in a weaker hue. 
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Figure 9 
Comparative effects of kerosene and gas oil on level of diffusion 


mm 
20, 


2 ee 
Contaminant {%) 


Figure 10 
Continued diffusion after 24 hours’ air-drying 


3. Comparative effects of kerosene and gas oil contamination. 

i) Level of diffusion. For gas oil the level of diffusion is approximately twice as 
high for the same amount of kerosene (Figure 9). 

ii) Hue. Owing to its slow evaporation, gas oil gives a ycllow hue, as against 
orange for kerosene. 

ili) Diffusion after completion of test. An additional reading taken after 24 hours’ 
air-drying showed continued diffusion of the dye. Figure 10 shows the 
changes in diffusion level for this period, after which the process is fi- 
nally arrested. It can be seen that this continued diffusion is smaller for 
pure gasoline but considerable for gas oil, thus providing an additional 
means of identifying the latter. 
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Figure 11 
Downward diffusion 
(1) Pure gasoline. (2) 0.5% gas oil. (3) 1% gas oil. (4) 2% gas oil. 


iv) Downward diffusion. This phenomenon was observed in kerosene and gas 
oil-contaminated gasoline, but not in pure gasoline: (Figure 11). 


4. Quantitative determination of the admixture. 


This is based on colorimetric comparison of the chromatograms of the tested 
sample and a standard solution of the same gasoline with given admixture of con- 
taminant under identical conditions of temperature and humidity. An example is 
shown in Figure 12. a 


SUMMARY 


Small amounts of keresone and gas oil in M.T. gasoline are detectable by one of 
the following methods: 


a) Semimicrodistillation. 
b) Paper chromatography. 
a) Semimicrodistillation. Preliminary distillation by the A.S.T.M. D-86 method 


up to 94%, and distillation of 5 ml of the residue in a 10 ml Wurtz flask, permits 
accentuation of the rise in F.B.P. due to the small admixture of contaminant. 
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Figure 12 
Identification of contaminant 


(1) Chromatograms 19, 32, 35(2), 42— pure gasolines (concentrated dye) with different F.B.P. 
(different levels of diffusion). (2) Chromatogram 30-1 “% kerosene contamination (uniform orange hue) 
(3) Chromatogram 35(1) — 1% gas oil contamination (uniform yellow hue). 


b) Paper chromatography. The chromatographic method permits: 

I. Qualitative detection of kerosene or gas oil admixtures of 0.5°%% and above, 
indicated by the higher level of diffusion and changed hue of the added dye. 

2. Identification of the contaminant according to the hue of the chromatogram. 

3. Quantitative determination of the admixture. 


A comparison of the standard methods and the two proposed by the authors is 
given in Figure 13. The proposed methods are especially important in view of the 
low F.B.P. of the currently marketed local M.T. gasoline (30°C below the 
standard specification), ruling out detection by the conventional standard methods. 
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Figure 13 
Comparison of methods in respect of accentuation of differences 
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BOOK REVIEWS 


MODERN MATHEMATICS FOR THE ENGINEER, edited by E. F. BECKENBACH, (Univer- 
sity of California Engineering Extension Series, xx+514 p. McGraw-Hill, New York, 1956), $7.50 


This book is best described in the words of two of the authors: “(This) series of chapters... is 
devoted to the contributions that mathematics can make nowadays to engineering” (D. W. Leh- 
mer, p. 480). 

“The chapters of this book fall roughly into three classes. Some deal with general theorems upon 
which statements concerning the general behaviour of physical systems can be based.... others 
deal with the formulation of problems to facilitate effective mathematical attacks. The rest.... deal 
mainly with questions of finding detailed numerical solutions for specific problems ....”’ (Ch. B. 
Tompkins, p. 448), 

The title of the book, however, is somewhat misleading. It is not a reference-book for engineers, 
but a stimulating exposition of applications of modern ideas in mathematics to a wide field of prob- 
lems—from ballistic tables to the Bush analyser and spurious business cycles. 

The book opens with a preface by the editor and an introduction by R. Weller, and consists 
of three parts: Mathematical Models, Probabilistic Problems and Computational Considerations, 
comprising 19 separate chapters of varying length, from 6 pages (No. 2) to 39 pages (No. 13 and 
16). The level is equally variable: some of the articles begin with definitions of fundamental con- 
cepts, while others plunge right into the latest advances. 


1, S. Lefschetz - Linear and Nonlinear Oscillations 


A very readable account of damped and forced oscillations, linear and nonlinear differential 
systems and their cycles. Duffing’s method of approximation is discussed in some detail. 


2. R. Bellman — Equilibrium Analysis: The Stability Theory of Poincaré and Liapunov 


A brief discussion of some stability problems related to linear differential equations, differential- 
‘difference equations and partial differential equations. 


3. J. W. Green - Exterior Ballistics 


The author describes the types of problems encountered in ballistics, indicates the coordinate 
system appropriate to each type, and mentions the special difficulties involved in increased ranges 
due in part to the Coriolis effect. He also discusses aerodynamic forces, points out the conside- 
rable amount of computation entailed and shows how ballistic and firing tables are constructed 
and the use of differential corrections. References are up to date. 


4. M.R. Hestenes - Elements of the Calculus of Variations 

A concise and clear introduction to the calculus of variations. Beginning with the brachisto- 
chrone and geodesics the author derives and proves in full the Euler conditions including, rather 
unusually, the Weierstrass and Legendre conditions. This is followed by some applications, a de- 
tailed discussion of the Hamiltonians, isoperimetric problems, etc. The chapter is concluded with 
a general formulation of the problem of Bolza and the necessary conditions for solution. 


5. R. Courant - Hyperbolic Partial Differential Equations and Applications 


With his usual clarity, the author discusses Brownian motion, why signal transmission is only 
possible in three dimensions, shock-waves and finite-difference methods. 


6. M. Schiffer - Boundary-value Problems in Elliptic Partial Differential Equations 


After showing that restriction to analytical functions only is impossible, Prof. Schiffer discus- 
ses properly posed boundary conditions for the Laplace equation and shows the connection bet- 
ween the Green, Neumann and Robin functions and the kernel function. This is followed by the 
maximum principle and variation formulae for Green’s functions. 


7. I. S. Sokolnikoff — The Elastostatic Boundary-value Problem 


The author describes Russian methods of solving two-dimensional elastotatic problems, ex- 
plains how a biharmonic boundary problem and a Dirichlet problem for cylinders are arrived 
at, and uses complex variables in solving them. 
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8. N. Wiener — The Theory of Prediction 
The author formulates and proves the mean and Birkhoff ergodic theorems, and solves the pre- 
diction problem. 


9. H. F. Bohneblust - The Theory of Games 
Zero-sum two-person games and the max-min = min-max relation for mixed strategies are 


discussed. . 
Examples include an application to economics and a poker-like game with possible chance moves. 


10. G. W. King - Applied Mathematics in Operations Research 

This makes easy reading. The author deals first with deterministic models (including a reference 
to linear programming, and the use of duals) and then with probabilistics models, including the 
newsboy’s inventory problem (with an interesting comparison of results arrived at by means of 
Poissons’ law, Monte-Carlo methods, and the introduction of a goodwill factor). This is followed 
by the “good customer” problem involving Laplace transforms, the use of the prediction theory 
and the Markoff processes. Generally speaking, this chapter comprises topics dealt with in some 
of the other chapters and illustrates the combined use of their conclusions. 


11. R. Bellman — The Theory of Dynamic Programming 

The term ‘‘dynamic”’ serves to emphasise the essential roles of time and of the sequence of deci- 
sions. The first problem is that of purchasing equipment, involving a given proportion of salvage. 

The author shows how recourse to classical calculus can be avoided by the use of infinite-stage 
approximation, proves existence and uniqueness theorems for the solution of maximum prob- 
lems in N-dimensional regions, and stresses the important fact that a solution either can be arrived 
at by guessing in function space or in policy space. This is followed by the gold-mining equation 
and a typical bottleneck problem (both in the form of continuous versions) and finally games of sur- 
vival. A comprehensive bibliography (55 titles) is given. 


12. G.W. Brown — Monte-Carlo Methods 


The author shows the use of sampling techniques for the evaluation of integrals, stressing its 
superiority over ordinary computation procedures for multiple integrals. This is followed by a dis- 
cussion of diffusion processes, and of sampling methods from specified distributions. 


13. L.A. Pipes — Matrices in Engineering 


The author discusses matrix algebra, calculus of matrices, integration of linear differential equati- 
ons and applications to elastic structures (including Castigliano’s theorem), and electrical engi- 
neering (especially three-phase networks and the impedance matrix) and vibration problems. A table 
is given of transmission matrices for fundamental four terminal structures (introduced by the author). 


14, J. L. Barnes — Functional Transformations for Engineering Design 
After considering the problems of design of an autopilot and its different stages, the author 
classifies them in three types: 
1) System-analysis (where driving functions boundary conditions and the system equations 
are known, and the responses are to be found). 


2) Instrumentation problem (where responses and system equations are known, the initial con- 
ditions are to be found. 


3) Design-problem (where both the initial date and responses are known, and the system equations 
are to be found). 


He gives examples of integrodifferential equations representing certain analogue systems, such as 


5 d (Li) 
v ()— —— —Ri—s | idt =0 
= — j 


and shows the various methods of transformations used. 
Two tables of integral transformations (Laplace, Fourier, Mellin, Hankel and Bessel) are included, 
Finally the author discusses the choice of the appropriate functional transformations. 
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15. E.F. Beckenbach — Conformal Mapping Methods 


A simple and clear presentation of the principles of equal-area, conformal, conical and other 
maps, including a discussion of stereographic projection conformality and the complex-variable 
theory, Laplace’s equation, the Riemann mapping theorem, the Osgood-Carathéodory theorem 
schlicht maps, Moebius transforms, and details of certain conformal mapping methods. 


16. Ch. B. Morrey Jr. — Nonlinear Methods 


_ The author describes Newton’s method for solving n equations in n unknowns and its modifica- 
tions, proves an inverse mapping theorem by vector methods, and discusses linear, normed-linar 
and complete spaces. He then deals with functionals and variational problems, and shows how to 
reduce non-linear systems to linear systems. Detailed examples of computatio nof numeticaw! 
lutions are given in an Appendix. The first example is the system 


4+ x+ y—x2+2xy+3y2 =0 
(a) 
14+ 2x—3y+x2+ xy—2y? =0 


solved by Newton’s method modified and then by the method of continuity, introducing the sys- 
tem 


(4+x+y+y(—x? + 2xy+ 3y) =0 


(p) 
(i+ 2x—3y)+ y@? + xy—2y?) =0 


which reduces to (a) for 2 = 1 and toa linear system for A = 0. The method consists in continuing 
the solutions for y = 0 by numerical integration to 2 = 1. 
Finally,” the author computes the solutions of 


1 
pecticn [y’(x)]2 dx = min 


oy 


10) 


given, y(0) = 0, y(1)=2/2using a modified Rayleigh-Ritz method. This involves a conside- 
rable amount of work but gives a very good approximation. 


17. G. E. Forsythe — What are Relaxation Methods? 

In the author’s own words “‘the question is an earnest one’”’—and so is the answer. The chapter 
gives the description of several methods, including an analysis of the cyclic relaxation method 
for a linear system of equations, a study of overrelaxation and its effect upon the speed of conver- 
gence, and an example of application to a non-linear problem. The bibliography lists 46 well- 


classified items. 


18. Ch. B. Tompkins — Methods of Steep Descent 

The author first describes the Kaczmarcz method for linear equations, then gives the differen- 
tial equations of steepest descent, and discusses the best possible lenght of each step, the Hestenes- 
Stiefel method and constraints. This is followed by extension to spaces of infinitely many dimen- 
sions and a study of the difficulties with regard to convergence of the solution. 

In conclusion, there are remarks on the close relationship between the method of steep descent 
and that used by the physical system itself in solving the problem. 


19. D.H. Lehmer — High-speed Computing Devices and their Applications 
Digital and analogue computers are compared, and sources of error discussed (truncation and 
round-off errors). This is followed by inversion of matrices and eigenvalue problems. 


ee 


The book concludes with name and subject indexes and pictures of the contributors are given 
on the back of its jacket. 
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As a final appreciation, it can be said that the editor has certainly made the right choice in every 
subject. Each aspect of applied mathematics is taken up by a recognised expert in that field and 
in most cases the level is equally high by both scientific and didactic criteria. 

This is indeed a book for every good research engineer, and should also be read by every mathe- 


matician. A. GUTWIRTH 


YEASTS. K. ArimA, W. J. NICKERSON, M. Pyke, H. SCHANDERL, A. S. SCHULTZ, A. C. THAYSEN 
AND R. S. W. THorRNE. Dr. W. Junk, Publishers, The Hague. 1957. 246 pp. $ 7.00. 

This is the first book of a series Biologia and Industria (edited by Roman of Adelaide, Genevois 
of Bordeaux, Parodi of Buenos Aires, Quastel of Montreal and Waldenstr6ém of Stockholm, which 
fills a lacuna, as no publication so far has tried, for biology with industry, while similar links to 
industry for chemistry and physics are available. The scope is international, and this aim has been 
well reached, as the present volume shows. The closer contact between scientist and technologist 
can only be to mutual benefit. This volume, Yeasts, illustrates well the difficulties involved in such 
a survey, which includes e.g. for the baker’s yeast the culture media, the biology, chemistry and 
physiology of that yeast and finally its industrial production. Similar is the arrangement for brewer’s 
yeast, for wine and fruit yeasts, for ‘sake’ (the eastern beer) and similar yeasts, for food- and fodder- 
yeasts, for medical, nutritional and technical yeast preparations. 

The book is well presented, as would not otherwise be expected. from the publisher. The scientific 
surveys are comprehensive. We find mentioned the names of a number of workers of The 
Hebrew University such as A. Fodor, J. Leibowitz and S. Hestrin. We feel it must be a real 
comfort for the yeast scientist to have this general survey which fully includes the technology of 
every branch, as it must be for the technologist to have this competent manual. We look forward 
impatiently to the further volumes of the series. 

F. S. BODENHEIMER 


INTRODUCTION TO CHEMICAL ENGINEERING THERMODYNAMICS, by J.M. SmitH 
and H.C. vAN Ngss, 2nd ed., 1959, McGraw-Hill Book Company, Inc., N. Y., 


Since S.L. Soo published his ‘“‘Thermodynamics of Engineering Science” in 1958, one might 
have hoped that the barriers between ‘‘Mechanica]’”’ and “‘Chemical” Engineering Thermodynamics 
would sooner or later break down. However, the present authors still believe, according to the 
preface, in ‘the justification of a separate text for chemical engineers’’. The plain fact is that at the 
utmost the last hundred out of a total of nearly 500 pages have a definite connection with chemistry, 
while the First and Second Laws, Fluid Properties and Flow Processes, Heat Engines and Refri- 
geration are the common property of all engineers (and others, too) even when a compression 
problem refers to methane instead of air or steam. 

With these reservations, the work is a highly competent textbook written in a clear and lucid style, 
with only occasional lapses into the colloquial, such as the ‘‘considerable investigation by famous 
workers” (p. 413) which led to the Third Law. The latter would have deserved more than a footnote 
(p. 195) which seems to attribute the authorship to Planck. 

Ten years after the 9th General Conference for Weights and Measurements recommended the 
gradual introduction of the Giorgi (MKSA) system, an up-to-date work should also have something 
to say about the units of measurements. At a time when every schoolboy of average ntelligence 
knows about the variability of the gravity acceleration in space flight (based on flow processes!), 
velocity formulae containing the gravity constant are far from convincing though conventionally 
correct. Here would have been the proper place to include references to the MKSA system and 
the ensuing elimination of g. 

There are several minor matters which careful editing ought to have eliminated: the relationship 
between adiabatic and isentropic processes is not clarified; the reference to the latter (p. 192) in 
the index is obviously misprinted. In the chapter on Refrigeration the authors rightly state that 
eq. (10-22) can be solved from a knowledge of the enthalpies, but not every reader will recall that 
an entropy chart for air is included a hundred pages before this equation. Similarly, the somewhat 
liberal statement that “the thermodynamic treatment of absorption refrigeration is not comli- 
cated” could be substantiated by the use of an enthalpy-concentration diagram (pp. 158-162) which 
should include the vapour region. 

All this, however, will not alter the fact that a student mastering this book will acquire an excel- 
lent knowledge of thermodynamics, even beyond the scope of chemistry. 

R. LANDSBERG 
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